Hadronic decays of mesons and baryons in the Dyson- Schwinger approach 
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We study hadronic decays of mesons and baryons in the context of the Dyson-Schwinger equa- 
tions of QCD. Starting from a well-established effective interaction in rainbow-ladder truncation, 
we consistently calculate all ingredients of the appropriate decay diagrams. The resulting strong 
couplings are presented as functions of the quark mass from the chiral limit up to the respective 
decay thresholds. In particular, we investigate the p — > nn and for the first time the A — > iV7r 
transitions. Both meson and baryon results compare well to available lattice QCD results as well 
as experimental data and present the first step towards a comprehensive covariant study of hadron 
resonances in the Dyson-Schwinger approach. 
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I. INTRODUCTION 

In hadron physics the strong interaction dominates 
the decay width of a resonance if appropriate hadronic 
channels are open. Thus, strong decays are an issue of 
paramount interest and any study of hadrons as mere 
bound states is necessarily incomplete if the states under 
investigation lie above the relevant thresholds. 

Consequently, strong decays of hadrons in terms of 
composites of quarks and gluons have been studied from 
the beginning, e.g., in quark models, and various de- 
cay mechanisms, formulations, and levels of sophisti- 
cation have been employed over the years. We ex- 
emplarily mention harmonic-oscillator models [1-3], the 
elementary-emission model [4-7], flux-tube breaking [8- 
10], the quark-pair creation or 3 Pq mechanism [11-13] 
and a simple coupled-channel formulation of meson res- 
onance properties in relativistic Hamiltonian dynamics 
[14-17]. These approaches range from the aforemen- 
tioned non- and semirelativistic as well as fully Poincare 
invariant quark-model calculations with various inter- 
actions to reductions of a Bethe-Salpeter treatment of 
hadrons, e.g., [18-21]. Ideally one would have a con- 
sistent and comprehensive coupled-channel calculation 
of meson and baryon resonances, see e.g. [22]. How- 
ever, to consistently implement an analogous approach 
within QCD is challenging. Recent lattice studies provide 
promising progress [23-29] ; nevertheless it is encouraging 
that reasonable results can already be obtained without 
working with resonances from the very beginning. 

In the Dyson-Schwinger-equation (DSE) approach to 
QCD two-body systems are described by the Bethe- 
Salpeter equation (BSE). Hadronic transition processes 
were considered in [30, 31], and an exploratory study 
of the strong decay of light vector mesons to two pseu- 
doscalars was performed in this framework some years 
ago [32] . Here we revisit this study and investigate in de- 
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tail both model-parameter and quark-mass dependence 
of the results, and we compare to experimental data and 
analogous calculations in lattice-regularized QCD. 

As the next important step towards a comprehensive 
investigation of the entire hadron spectrum, we gen- 
eralize this calculation to a strong baryon decay and 
calculate the hadronic coupling of the A— baryon to 
Nir. The A(1232) resonance plays an important role 
in pion-nucleon scattering and pion photoproduction ex- 
periments; see [33, 34] for recent reviews. Its experi- 
mental width mainly owes to the decay into a nucleon 
and a pion, whereas the electromagnetic A —> Nj de- 
cay channel is considerably suppressed. The correspond- 
ing strong coupling constant g ANV has been studied in 
various approaches, e.g., in the quark model [35, 36], 
in meson- exchange models [37, 38], via light-cone sum 
rules [39, 40], and recently also in lattice QCD [41]. 
The A— baryon is an essential component in a realis- 
tic description of meson-baryon interactions in nuclear 
physics, for instance via chiral effective field theories, and 
the computation of the AiV7r vertex from its underlying 
quark-gluon dynamics is an important task in hadron 
physics. 

The present study is motivated by recent successes in 
the implementation of the Dyson-Schwinger approach to 
various aspects of hadron phenomenology. In particular 
we make use of a well-established effective model setup in 
rainbow-ladder (RL) truncation, together with the con- 
sistent extension to a quark-diquark picture for baryons. 
Successful applications of this setup pertain to various 
observables: meson spectra and leptonic decay constants 
of pseudoscalar and vector mesons were studied over the 
full range of quark masses from chiral limit to bottomo- 
nium [42, 43], and recently the feasibility of such meson 
studies for any spin has been demonstrated [44] . Electro- 
magnetic properties of pseudoscalar and vector mesons 
involve a consistent construction of the electromagnetic 
interaction process via triangle diagrams analogous to the 
one used here [45-48] . They have proven to be an excel- 
lent example for the importance of correctly implement- 
ing the symmetry properties of the underlying theory in 
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numerical calculations. 

The extension of the approach to baryons can be sim- 
plified via the introduction of diquarks [49-51]. The cor- 
responding investigations of baryons in a covariant quark- 
diquark setup have already undergone considerable devel- 
opment [52-58] and very recently culminated in the first 
genuine three-quark treatment of nucleon and A masses 
and nucleon electromagnetic form factors in the covariant 
Faddeev framework [59-62]. 

This article is organized as follows: in Section II we 
collect the necessary building blocks of the approach; the 
p — > tttt and A — > Nt: transition diagrams are worked 
out in Section III; Section IV contains our results for both 
meson and baryon sectors; and we conclude in Section V. 
Decay-width formulas and the color-flavor traces of the 
decay diagrams are collected in the appendices. All cal- 
culations are performed in Euclidean momentum space 
and in the isosymmetric limit in Landau gauge QCD. 

II. BUILDING BLOCKS 

Our investigation of hadronic decays involves numeri- 
cal solutions of several integral equations whose proper- 
ties and solutions have been studied elsewhere. In partic- 
ular, we are concerned with the quark DSE, the meson 
and diquark BSEs and, in the context of baryons, the 
quark-diquark BSE. In the following we briefly review 
these equations and the properties of their solutions, and 
for each case we refer the reader to more detailed discus- 
sions in the literature. 



A. Truncation and effective interaction 

Numerical model studies of hadrons such as the one 
presented here necessitate a truncation of the infinite 
tower of DSEs. In the following we will restrict our- 
selves to the RL truncation which substitutes the fully 
dressed quark-gluon vertex with a bare vertex. Its coun- 
terpart in a hadronic bound-state equation is a gluon 
ladder kernel, i.e. a dressed iterated gluon exchange be- 
tween two quarks. The combined strength of the gluon 
propagator and quark-gluon vertex is then modeled by 
an ansatz. Phenomenologically important directions of 
improvement beyond RL involve the implementation of 
pseudoscalar meson-cloud effects but also other struc- 
tures in the quark-gluon vertex and the qq and qq kernels, 
see [63-65] and references therein. 

The RL truncation offers many advantages. It is sim- 
ple to implement, but at the same time allows for so- 
phisticated model approaches to QCD within the DSE- 
BSE context since it satisfies the axial-vector and vector 
Ward-Takahashi identities (see e.g. [45, 46, 66-71]). The 
axial-vector Ward-Takahashi identity is essential for the 
correct realization of chiral symmetry and its dynamical 
breaking in any model calculation. In particular, it im- 
poses constraints upon the construction of the integral- 



equation kernels. As the most prominent result, Gold- 
stone's theorem is satisfied [70] and one obtains a general- 
ized Gell-Mann-Oakes-Renner relation valid for all pseu- 
doscalar mesons and all current-quark masses [72, 73]. 
The vector Ward-Takahashi identity on the other hand 
is the guiding principle for the construction of consistent 
electromagnetic currents. 

In RL truncation the equations as presented in the fol- 
lowing subsections contain an essential model ingredient, 
namely an effective interaction Q which we choose from 
Ref. [74] as 

g(fc 2 ) Dnk* fcW 2n lm (1 - e- fc2 /^) 

4*Z* ' ufi + ln[r+(l + P/A| CD )2]' 1 > 

where k is the gluon momentum and Zi the quark renor- 
malization constant. This particular form has been em- 
ployed in many of the works listed as successes of the 
approach in the introduction. It provides the correct 
amount of dynamical chiral symmetry breaking as well as 
quark confinement via the absence of a Lehmann repre- 
sentation for the dressed quark propagator. Furthermore, 
it has the correct perturbative limit, i.e. it preserves 
the one-loop renormalization group behavior of QCD 
for solutions of the quark DSE. Following [74], we have 
A t = 1 GeV, r = c 2 — 1, Nf = 4, Aq^ = 0.234 GeV, and 
7 m = 12/(33 — 2Nf). The main motivation of this func- 
tion, which mimics the behavior of the product of quark- 
gluon vertex and gluon propagator, is of phenomenologi- 
cal origin. While currently debated on principle grounds 
(e.g. [75, 76]) the impact of its particular form in the far 
IR on meson masses is expected to be small (see also [77] 
for an exploratory study in this direction). 

The phenomenologically important regime, in partic- 
ular with respect to light meson (e.g., pion) properties, 
is the intermediate-momentum region, modeled by the 
Gaussian term in Eq. (1). D and co, in principle free 
parameters of the model interaction, can be used to in- 
vestigate certain aspects of both the interaction and the 
bound states in the BSE. In particular one can interpret 
D as an overall strength and co as an inverse effective 
range of the interaction. In the range to € [0.3, 0.5] GeV, 
the prescription D x lo —: Aj R = const, follows from 
fitting of the model parameters to ground-state proper- 
ties [74] and defines a one-parameter model, character- 
ized by a fixed infrared scale Air = 0.72 GeV. Appar- 
ently [43] an insensitivity of an observable to uo in this 
prescription is characteristic for a ground state, while or- 
bitally or radially excited-state properties show consid- 
erable dependencies on uo. In the present work, we also 
consider such a possibility and plot bands to indicate the 
model-parameter dependence of our results. 

B. Quarks 

The Dyson-Schwinger equation for the quark propaga- 
tor in rainbow truncation is illustrated in the left panel 
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FIG. 1. (Color online) Left panel: quark DSE (2) in rainbow truncation. Right panel: schematics of the meson BSE (3) and 
diquark BSE (5) in ladder truncation. 



of Fig. 1 and reads 

sip)- 1 = z 2 (ip+m ) + g{k 2 ) ^ y S(q) Y , (2) 

where S(p) is the renormalized dressed quark propaga- 
tor, p and k — p — q are the quark and gluon momenta, 
and J = J <i 4 g/(27r) 4 represents a four-momentum inte- 
gration. T£7& 2 with T% v = 5^ - k^k v /k 2 is the free 
gluon propagator, 7^ is the bare quark-gluon vertex, and 
G(k 2 ) is the effective interaction defined above in Eq. (1). 
Dirac and flavor indices have been omitted for simplic- 
ity and the factor 4 / 3 comes from the color trace. The 
bare current-quark mass mo is the input of the equation. 
The solution of Eq. (2) requires a renormalization proce- 
dure, the details of which can be found together with the 
general structure of the quark DSE in e.g., [72, 74]. 

S(p) is an important ingredient in all of the following. 
We note that its solution as an input for the various BSEs 
described below must be known in a parabola-shaped re- 
gion of the complex p 2 plane, the size of which is propor- 
tional to the mass of the respective bound state that it 
helps to constitute. As a result, a sophisticated numeri- 
cal approach is needed and we refer the reader to [78] for 
a description of our particular solution method. 

C. Mesons 

In our framework mesons with total qq momentum P 
and relative qq momentum p are studied via the meson 
BSE. Its structure in RL truncation is shown in the right 
panel of Fig. 1 and reads 

r M (p, p ) = -\J S(k 2 ) 5Vx M (g, P) i v , (3) 

where Tyi{k;P) is the Bethe-Salpeter amplitude and 
Xm{<1\P) = S(q + )TM{q; P) S(q-) is referred to as the 
Bethe-Salpeter wave function. The quark and anti- 
quark propagators depend on the (anti)quark momenta 
g+ = q + r)P and g_ = q — (1 — r/)P, where r\ £ [0, 1] is 
a momentum partitioning parameter usually set to 1/2 
for systems of equal-mass constituents (which we do as 
well). The quark-antiquark interaction kernel is given 



by a ladder dressed-gluon exchange, whose dependence 
on the gluon momentum is characterized by the same ef- 
fective interaction, Eq. (1), as in the quark DSE Eq. (2). 
The combined set of truncated equations (2) and (3) thus 
already by construction satisfies the axial-vector Ward- 
Takahashi identity. 

The general dependence of the meson amplitude on 
the various four-momenta can be written in terms of N 
covariant structures Tj reflecting the quark and antiquark 
spins, together with scalar components Fi, i = 1, . . . , N 
as 

JV 

r M (q-P)=Y,T l (P;q;i)F l (q 2 ,q-P,P 2 ) , (4) 
i=l 

where semicolons separate four-vector arguments, and 
N = 4 for mesons with total spin J — and N = 8 
otherwise (see e.g., [44]). A detailed account of the Tj for 
pseudoscalar and vector mesons as needed in our calcu- 
lation can be found in Ref. [43]. 

The components are scalar functions of their three 
scalar arguments: the total momentum squared P 2 , the 
relative momentum squared q 2 , and an angular variable 
q ■ P. Note that for an on-shell amplitude P 2 = —M 2 
is fixed, while one artificially varies P 2 in the solu- 
tion process of the homogeneous BSE (see e.g., [79, 80], 
where also all necessary details on the numerical solution 
method can be found). In the corresponding inhomoge- 
neous vertex BSE on the other hand, one would have P 
and therefore also P 2 as a completely independent vari- 
able (see, e.g. [79, 81, 82]). Thus, the on-shell scalar 
components Fi(q 2 ,q-P, P 2 ) effectively depend on the two 
variables q 2 and q ■ P. The latter can be parameterized 
by the variable z e [—1,1] related to the cosine defining 
the angle between P and q. With such a reparameteriza- 
tion in mind, the components Fi can be expanded further 
in Chebyshev polynomials, which leaves Chebyshev mo- 
ments of the Fi as sole functions of q 2 (for details and an 
illustration of Chebyshev moments, see [72, 83]). 

For our case very few Chebyshev moments are suffi- 
cient to produce converged results. However, in the con- 
text of the decay processes we note that, as a result of 
the kinematics in the triangle diagram as shown below, 
one needs to know the Bethe-Salpeter amplitudes in a 
certain region for the relative momentum squared p 2 in 
the complex plane. We achieve this by a continuation of 
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the relevant Chebyshev polynomials into the complex p 2 
plane via a Taylor-expansion technique up to 4th order 
which yields a converged result. 

D. Diquarks 

The relevance of diquark degrees of freedom in hadron 
physics has been reviewed in [49], and diquarks are in 
many respects conceptually similar to mesons. In our 
setup, diquark correlations appear as structures in a 
quark-quark system whose properties can depend on the 
truncation or the effective interaction. In particular, in 
RL truncation diquarks appear as timelike poles in the 
quark-quark T— Matrix, which is an unphysical result 
since diquarks are not color singlets but elements of an 
antisymmetric color antitriplet. This has been identified 
as a truncation artifact, i.e., diquarks disappear from the 
physical spectrum beyond RL truncation [84] . Neverthe- 
less, the significance of diquark correlations as binding 
structures within baryons has become apparent in vari- 
ous baryon form- factor studies, see e.g. [56, 58, 62, 85]. 
Moreover, the diquark concept has received support from 
investigations of diquark confinement in Coulomb-gauge 
QCD [86]. 

In our particular case the diquark BSE in RL trunca- 
tion reads 

T D (p 7 P)C = -* J g(k 2 ) T f^ X D(q,P)C 1 ^ (5) 
1 

where C is the charge-conjugation matrix, Yjj is the di- 
quark amplitude, xd(<Z! P) C — S(q+.)Tn(q; P) C S(q-), 
and for all practical purposes the only difference to the 
meson BSE Eq. (3) is the color factor. 

In the description of baryons as quark-diquark bound 
states and also for the description of baryonic transitions 
as given below one needs to know not only the diquark 
amplitudes but also the diquark propagator. A defining 
equation for this propagator can be consistently derived 
from the two-quark Dyson equation and reads schemati- 
cally [59] 

D- 1 = trj T D ST D S T -tv J I TdK^Td, (6) 

which still contains on-shell amplitudes T d resulting from 
Eq. (5). To obtain the propagator D for general P 2 ap- 
propriate ansatze for the off-shell amplitudes are cho- 
sen [59]. K is the same quark-quark interaction kernel 
that enters Eqs. (3-5), and a bar on an amplitude denotes 
charge conjugation: T(q,P) = CT(—q,—P) T C T . 

E. Baryons 

In this work baryons are interpreted as bound states 
of a quark and a diquark, which reduces the three-quark 
problem to an effective two-body problem derived via 




FIG. 2. (Color online) The quark-diquark BSE, Eq. (7). 



omission of three-quark interactions and a pole ansatz in 
the quark-quark T-matrix, see e.g., [52-54, 59]. The in- 
teraction in the resulting equation is an iterative quark 
exchange, where in every iteration step the spectator 
quark and one quark inside the diquark exchange roles. 
This quark-diquark BSE is illustrated in Fig. 2 and reads 

r%(p,P)=cWj K^ DQ S(k q )D^'(k d )T^(k,P), (7) 

k 

where Tg are the quark-diquark amplitudes of the respec- 
tive baryon. Their general Dirac-Lorentz structure is de- 
scribed e.g. in Refs. [52, 59] and their decomposition in 
terms of covariant basis elements and Lorentz-invariant 
components is analogous to the previously discussed case 
of a meson amplitude, i.e., Eq. (4) and below. The quark- 
diquark exchange kernel is given by 

^Q-DQ = rg(*V, k d ) S T (q) TUPnPd) , (8) 

with p qJ i and fc 9j£ ; being the external and internal quark 
and diquark momenta and p r , k r the relative momenta 
that enter the diquark amplitudes, cf. Fig. 2. 

The treatment of light baryons such as the nucleon 
and A in the quark-diquark approach usually retains the 
lightest diquark degrees of freedom, i.e., scalar and axial- 
vector diquarks. The A baryon then involves only axial- 
vector diquark correlations whereas the nucleon contains 
both. As a consequence, the ANir coupling will in- 
volve axial-axial contributions as well as axial-scalar di- 
quark transitions. The superscripts in Eqs. (7-8) ac- 
count for both kinds of diquarks and an implicit sum 
over these indices is understood: I Q are scalar (a = 0) 
or axial-vector diquark amplitudes (a = 1 ... 4) obtained 
from their respective diquark BSEs (5), and D 00 and 
denote the scalar and axial-vector diquark prop- 
agators, respectively. The color-flavor trace in Eq. (7) 
reads = -c^ = -1/2 and C (° Q ) = C ( Q °) = Vs/ 2 

for the nucleon, whereas in the case of the A it is given 
by = — 1. For details on the solution of the quark- 
diquark BSE we refer the reader to Refs. [59, 87]. 

III. HADRONIC DECAYS 

Given a certain truncation of the DSE-BSE system, 
one can compute observables involving various currents 
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FIG. 3. (Color online) The pirn triangle diagram, Eq. (10). 



through a consistent construction of the relevant invari- 
ant transition matrix elements. In our case we consider 
transitions between three hadrons and need a quark-level 
picture of the corresponding matrix elements. In the me- 
son case, for RL truncation one arrives at a so-called tri- 
angle diagram [32], depicted in Fig. 3, which corresponds 
to a generalized impulse approximation. Analogous di- 
agrams are used at this level for, e.g., meson electro- 
magnetic form factors [45]. For baryons an appropriate 
construction is also possible and given below. 



More generally, the pnir transition matrix element cor- 
responds to the coupling of the pion to an external 
vector current. Specifically, if the p— meson amplitude 
that appears in Eq. (10) were replaced by the dressed 
quark-photon vertex, evaluated at arbitrary momentum 
squared Q 2 , the respective triangle diagram would con- 
stitute the pion's electromagnetic form factor: A^^ = 
2P fJ 'F 7r (Q 2 ). As the photon fluctuates into p°, the 
dressed quark-photon vertex, obtained from its inhomo- 
geneous BSE, self-consistently develops a p— meson pole 
whose on-shell residue is proportional to the p— meson 
amplitude [45]. Such a purely transverse compo- 
nent is an important ingredient in various hadronic form- 
factor studies where it contributes typically ~ 50% to n, 
N and A squared electromagnetic radii [45, 48, 85, 88]. 
Consequently, the residue of the pion form factor at the 
meson pole Q 2 



P 



Fir(Q = -m 2 ) 



m p is proportional to g p7r7r : 

fp Wlp <?p7T7r 



Q 2 



V2 



(11) 



In the present study we are primarily interested in 
the value of g p7T7: on the mass shell Q 2 = — m 2 . In 
a covariant formalism such as ours, one may simply 
choose a frame of reference — in our case the rest frame of 
the decaying particle — by setting Q = (0, 0, 0, im p ) and 
P = (0, 0, K, 0), where n 2 — m 2 /4 — m 2 . Then, together 
with Eq. (10) , Qp-ir-n is straightforward to evaluate numer- 
ically. 



A. Mesons: p — > irn 

For the meson sector we investigate the p — > irir tran- 
sition which, in terms of Lorentz quantum numbers, cor- 
responds to an VPP-vertex. As in the usual kinematical 
setup for a two-body decay one has the total p— meson 
momentum Q and the relative momentum of the pion 
decay products P = (P2 — Pi)/2. All mesons are on- 
shell, i.e. Q 2 = —m 2 and P 2 = P 2 = —m 2 , which entails 
P ■ Q = 0. Due to the transversality of the p-meson the 
most general Dirac-Lorentz structure of the transition 
can be parametrized as 



A p 

pTTTT 



2P m 



9p7T7T ■ 



(9) 



where g p ^-n is its dimensionless coupling constant. The 
corresponding triangle diagram is illustrated in Fig. 3 and 
reads 

A^ = tr f r«(q 2 ,P 2 )S(q + P)Ti r (q 1 ,P 1 )xZ(q,Q), 

(10) 

where the are the (charge-conjugated) on-shell pion 
amplitudes, i.e. the canonically normalized solutions of 
their homogeneous BSEs, S is the renormalized dressed 
quark propagator obtained from its DSE, and Xp is 
the p— meson wave function defined in connection with 
Eq. (3). The traces in color and flavor space yield a fac- 
tor 3a/2 in front of the integral, cf. App. B. 



B. Baryons: A — > Ntt 

In the baryon case, the coupling of an on-shell nucleon 
and A— baryon to a pseudoscalar current is described by 
the pseudoscalar transition form factor Gantt(Q 2 )- The 
A-baryon is a spin-3/2 particle and is thus represented 
by a Rarita-Schwinger spinor. Combined with the re- 
striction to positive energies for the nucleon and A, the 
most general Dirac-Lorentz structure of the respective 
interaction vertex is given by 



A M 

A 7V7T 



G AN .(Q 2 ) ^- A+(P/) P^(P,) , (12) 

Z1V1-- 



N 



where Pi, P/ are the incoming A and outgoing nucleon 
momenta, Q = Pf — Pi is the pion momentum, and 

Ag(k) = (1 + #)/2 is the positive-energy projector of 
the baryon with k being the respective normalized mo- 
mentum Pi or Pf. The Rarita-Schwinger projector of the 
A reads 

TO=A+(i)fe"-i^j, (13) 

where is a transverse projector with respect to k and 
the 7— matrices 7^ are transverse to k as well. 

At the quark level, the coupling of nucleon and A 
to a pseudoscalar current is represented by the quark- 
pseudoscalar vertex T 5 . The latter satisfies an inhomoge- 
neous BSE which has the same structure as Eq. (3) except 






FIG. 4. (Color online) Decomposition of the ANn transition matrix element in the quark-diquark model, Eqs. (16-17). Seagull 
terms are neglected. 'SC and 'AX' denote the types of diquarks (i.e., scalar and axial-vector) that can appear in combination 
with the nucleon and A bound-state amplitudes. 



for an additional inhomogeneous term on the right-hand 
side [72]: 

r 5 (p, Q) = z i i l5 ~^J g(k 2 ) TL^ X5 ^ Q) i v ■ (14) 



Its residue at the pion pole Q 2 = —m 2 is proportional to 
the homogeneous pion amplitude r„. (cf. also [72, 82]): 



r 5 (p,Q) 



Z 2 2m Q 2 



1 



r„(p,Q), (15) 



where f„ is the (calculated) pion decay constant and mo 
is the bare current-quark mass that enters the quark 
DSE (2). The solution for T 5 provides an off-shell ex- 
pression for the pion amplitude rV and thereby allows 
to compute the pseudoscalar transition form factor at 
spacelike momenta Q 2 > 0. Then, Gan-k(Q 2 ) corre- 
sponds to the form factor obtained from T„ (p, Q) in 
Eq. (15), where the pion pole as well as its residue are 
removed, and its on-shell value is the AiV7r coupling con- 
stant GANiv(-ml) = gANv 

For the A — > Ntt system the construction analogous 
to Eq. (10) is more complex since there are diquarks as 
well as quarks present as constituents of the states in- 
volved in the transition. Moreover, the pion will not only 
interact with the quarks and diquarks directly but can 
couple to the quark-diquark kernel as well, i.e. impulse- 
approximation diagrams alone are no longer sufficient in 
studying the A — > Ntt transition. 

A systematic procedure to derive the coupling of a 
baryon to an external current is the gauging-of-equations 
method of Refs. [89-91]. In the context of baryon elec- 
tromagnetic form factors it has been applied to the 
quark-diquark model [92] as well as the three-quark ap- 
proach [62] . The starting point is to identify the current 
with the residue of the 'gauged' quark-diquark (or three- 
quark) T— matrix on the baryon's mass shell. Upon ex- 
ploiting the relation between the T—matrix and the ker- 
nel of the respective bound-state equation, the hadronic 
matrix elements of the current are obtained as a sum 
of diagrams that describe the coupling of the current 
to all ingredients at the constituent level, i.e. in our 



case to the quark and diquark propagators as well as the 
quark-diquark kernel. The procedure can be applied for 
mesons as well where, in the case of a rainbow-ladder 
quark-antiquark kernel, the triangle diagram of Fig. 3 is 
recovered. 

The generalization of the method from an electromag- 
netic to a pseudoscalar current, as well as different kinds 
of baryons in the initial and final state, is straightfor- 
ward. The resulting diagrams are displayed in Fig. 4 and 
involve impulse-approximation couplings to the quarks 
and diquarks as well as a coupling to the exchanged quark 
that appears in the quark-diquark kernel. In principle 
there would be further diagrams that contain seagulls, 
i.e., pseudoscalar couplings to the diquark amplitudes. 
For electromagnetic form factors such seagull contribu- 
tions are typically small [56, 59] but necessary to ensure 
electromagnetic gauge invariance; however, in the present 
case we do not consider them further. The ANtt transi- 
tion matrix element is then decomposed as 



A M 

A 7V7T 



A£ + A£ Q + A£ X , 
where the three contributions are given by 



(16) 



DQ 



A£x = 



r a v r 

1 TV A7T 1 



1 N d 1 



D a0 , 

D aa 'r?!f dp' . 



(17) 



s v^xl rg' s v p A D aa ' ' d^' . 



Here we suppressed the explicit momentum dependen- 
cies for brevity. The kinematics are analogous to the 
electromagnetic form- factor case and are described, e.g., 
in App. C.l of Ref. [58]. Y a N and T^ 1 are the quark- 
diquark amplitudes for the nucleon and A— baryon; Xtt = 
ST^S is the pion (off-shell) Bethe-Salpeter wave func- 
tion obtained through the pseudoscalar vertex (15); D afi 
is the diquark propagator; is the diquark ampli- 
tude; and is the vertex that describes the cou- 
pling of the pion to a diquark propagator. If an axial- 
vector diquark is involved, a, ft = 1 ... 4 are Lorentz in- 
dices. Scalar-axialvector transitions, originating from the 
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0.25 



mi [GeV 2 ] m£ [GeV 2 ] 

FIG. 5. (Color online) Left panel: Evolution of the pirn coupling with the pion mass squared. The experimental point is 
indicated by the star and the symbols denote lattice data from Refs. [23-28]. For better readability the point of the Gockeler 
group is shifted slightly to the right. Right panel: The decay width of the p-meson versus the pion mass squared. The width 
of the bands illustrates the dependence on lu (see text). 



scalar-diquark component in the nucleon, can only occur 
in A^Q and A^ x ; in that case: a = 0, and D 00 de- 
notes the scalar diquark propagator and T^L the scalar- 
axialvector transition vertex induced by the pion. 

We note that all ingredients of Eq. (17) are deter- 
mined self-consistently: Eqs. (2), (5-7) and (14) provide 
the dressed quark propagator, the scalar and axial- vector 
diquark amplitudes and propagators, the baryon ampli- 
tudes and the pseudoscalar vertex. The diquark-pion ver- 
tices are obtained in analogy to Eq. (10), i.e. through 
respective triangle diagrams. The color and flavor traces 
in Eq. (17) are worked out explicitly in App. B. 



IV. RESULTS AND DISCUSSION 



The building blocks of the p — > ttit and A — > Nir tran- 
sition matrix elements have been determined in the pre- 
vious sections and we proceed with computing Eqs. (10) 
and (16) numerically. The Lorentz-invariant coupling 
constants <7 p7r7r and g &N7T are extracted via appropriate 
momentum contractions and Dirac traces, cf. Eq. (A12). 
Since within our chosen truncations the ingredients of 
the equations for g p7nr and g& NT; are computed selfcon- 
sistcntly, the effective quark-gluon coupling defined in 
Eq. (1) is the only model parametrization with impact 
upon the results. We explore the sensitivity to this ansatz 
by varying the lo parameter from its central value lo = 0.4 
GeV which is indicated by the colored bands in the plots 
of this section. 



A. Mesons 

Our result for g p7V ^ at the physical u/<i-quark mass 
is shown in Table I and compared to the experimental 
point given by the PDG [93]. In Fig. 5 we plot <7 p7r7r as a 
function of the pion mass squared and compare to recent 
lattice results [23-28]. 

The following observations are important: first, both 
the magnitude and dependence of our results are in 
agreement with lattice results, including a slight under- 
estimation of the experimental value by about 15%. Sec- 
ond, the lo dependence of the result is small, which re- 
sults from the ground-state characteristics of the n- and 
p-meson amplitudes regarding their dependence on the 
relative-momentum squared. This second point also so- 
lidifies the result of [32] where only the central value of 
our band, lj = 0.4 GeV was used. 

Next, the very weak dependence on compares 
well to analyses using chiral perturbation theory [94- 
96], where one concludes from the strong phase-space 
dependence of the decay width that the coupling's mj- 
dependence should indeed be close to zero. Consequently, 
when plotting the p — ¥ rnr decay width in our approach 
as a function of m^, Fig. 5, the falloff can be almost com- 
pletely attributed to the phase-space factor in Eq. (All) 
and the width vanishes when the decay channel closes. 

While comparison to experiment is favorable, the small 
difference between our result and the experimental value 
for g p7TTT warrants some discussion. The present model 
calculation uses a simple truncation, but an effective in- 
teraction which is fitted to the pion mass and underes- 
timates the p-meson mass by about 5%. While the re- 
sulting kinematical mismatch is responsible for a small 
part of the difference, the main contributions are oth- 
ers. Even though one has to expect some effect from 




FIG. 6. (Color online) Left panel: the transition form factor Gantt as a function of the squared pion momentum Q 2 . The 
upper band corresponds to the full result and the lower band shows the impulse-approximation value where the pion couples 
to the quark line only. The width of the bands again shows the dependence on uj. The solid lines are the results of the 
calculation in the kinematically allowed regions whereas the dashed lines are the respective dipole fits, cf. Eq. (19). The pion 
mass shell Q 2 = —m 2 , is indicated by the vertical line and the star shows the experimental value gANn = Gantt{— m 2 ^). Right 
panel: Current-mass evolution of gANn (upper band) and Gantt(Q 2 = 0) (lower band). The vertical shaded area depicts the 
oj-dependent location of the threshold ALa = Mn + m-* . 



non-resonant corrections to RL truncation, the main im- 
provement would be a self-consistent treatment of the 
p meson as a resonance, i.e., an inclusion of an explicit 
7T7T decay channel in the interaction kernel of the vector- 
meson BSE, in which, of course, also a different (refitted) 
effective interaction would have to be used. However, 
such an approach is much more involved than the present 
one and clearly beyond the scope of this study, in par- 
ticular for the baryon case. In addition, the reasonably 
small difference to the experimental value of g p7r n gives 
reason to expect that for this particular transition the 
present approach is at least a reliable gauge for future 
studies and results. 



B. Baryons 

The decay width of the A— baryon is governed almost 
exclusively by the strong interaction, namely via the de- 
cay into the nucleon and a pion. The only other decay 
channel, the electromagnetic A — » Nj transition, has 
a branching fraction of less than 1%. Experimentally, 
Ta = 118(2) MeV [93], from which the corresponding 
coupling strength g &N7T = 29.4(2) can be inferred via 
Eq. (All). Different conventions that are commonly em- 
ployed in the literature read 



2M N 







m p 


Mjv 


Ma 




f/AJV7T 


This work 


0.14 


0.74 


0.94 


1.28 


5.20 


28.1 


Experiment 


0.14 


0.77 


0.94 


1.23 


5.98 


29.4 



TABLE I. Comparison of our summarized mass and coupling 
values at the physical point for the central value of the uj 
band to corresponding experimental data. The pion mass 
is fitted to experiment, all other numbers are predictions of 
the model with no further parameters adjusted or introduced. 
The masses are given in GeV; the coupling constants are di- 
mensionless. 



with g'^ a 2.16 and g'l Nn a 15.7 GeV- 1 . 

For general off-shell momenta the A — > Nir coupling 
is described by the pseudoscalar transition form factor 
G A n7t(Q 2 ) which we obtain from Eqs. (16) and (A12). To 
compute the Q 2 — dependence of the form factor we work 
in the Breit frame where the pion momentum is given by 
Q = (0,0, |Q|,0). This has the advantage that the rela- 
tive momenta in the baryon amplitudes are real and no 
continuation into the complex plane is necessary. How- 
ever, due to the difference in the nucleon and A masses, 
the singularity structure in the quark and diquark propa- 
gators imposes kinematical constraints on the accessible 
Q 2 region from both below and above. Hence, to obtain 
the form factor g &N7r — G AJV7r ((3 2 = — m 2 .) at the pion 
mass, we fit our results at spacelike Q 2 with a dipole 
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form: 



GDipolc(Q 2 



G A JV7r(0) 

(1 + Q7AJ) 2 



(19) 



where G AN 7r(0) and are free fit parameters. 

Our result for the transition form factor G &N7T (Q 2 ) at 
the physical u/d mass is shown in the left panel of Fig. 6. 
Its computed value in the kinematically allowed range is 
plotted as a band with solid margins, where the width of 
the band corresponds to the model parameter u, whereas 
the fit results in the unaccessible region are shown as 
dashed lines. The resulting value of the strong coupling 
constant, g AN7T — 28.1, is remarkably close to the exper- 
imental number. We also note that the Q 2 — evolution of 
G AN7T (Q 2 ) is in good agreement with the lattice data of 
Ref. [29]. 

In general the dipole fits work very well and pro- 
vide an adequate representation of the form factor 
Gantt at spacelike values of the squared pion momen- 
tum. As described in connection with Eqs. (14-15), the 
Q 2 — evolution of the form factor is governed by the pseu- 
doscalar vertex T§ which includes all pseudoscalar-meson 
poles, i.e. both the pion's ground state as well as its exci- 
tations. While the pion ground-state pole and its residue 
were removed from Eq. (15) to obtain G AJV7r , the remain- 
ing excited states are still encoded in the vertex, hence 
the form factor G AJVlr must diverge at the respective pole 
locations. Indeed we find that the dipole mass in 
Eq. (19) roughly coincides with the mass of the first ex- 
cited state in the h channel which in RL trunctation, 
at the u/d mass and for the central to value, is obtained 
as m^* = 1.1 GeV [97]. 

The left panel of Fig. 6 also includes the (quark-) 
impulse-approximation contribution to G AJV7 r, i-e. the 
first diagram in Fig. 4 corresponding to Aq of Eq. (16). 
Here the diquark is merely a spectator and, since the 
A— baryon only involves axial-vector diquark degrees of 
freedom, only an axial-vector diquark propagator can ap- 
pear in that diagram. Fig. 6 shows that such a direct 
coupling to the quark provides roughly one third of the 
value of G AJV7r - The axial- axial contributions stemming 
from the second and third diagrams are comparatively 
small and contribute ~ 10% to the full result. The re- 
mainder owes in equal parts to the axial-scalar transi- 
tions that are generated from the transition vertex 
in the second diagram and the axial-scalar contribution 
r^Xw Tpj in the exchange diagram, cf. Eq. (17). 

Once again, the current-quark mass dependence of the 
transition form factor can be studied by varying the cur- 
rent mass in the quark DSE. That change will be re- 
flected in all ingredients that enter the transition ma- 
trix element. Similarly to the meson case, the over- 
whelming contribution to the mass dependence of the 
decay width r AJV7r comes from the phase-space factor 
in Eq. (All). This is especially conspicuous in the 
form factor G &N7r (Q 2 = 0) at vanishing pion momen- 
tum, shown in the right panel of Fig. 6, which is prac- 
tically independent of the current-quark mass. Similar 



features have been reported for TV and A electromag- 
netic form factors [58, 62]. The observation stays true 
for the Q 2 — evolution, i.e. G &N7r (Q 2 ) as well as its in- 
dividual contributions retain their shape throughout the 
current-mass range if they are plotted over a dimension- 
less variable such as Q 2 /Mfj. Considering Eq. (19), this 
means that the dipole fit works also well for higher quark 
masses since the mass of the excited pion also varies with 
the current-quark mass in a similar fashion as Mjv and 

On the other hand, the value of g AJV 7r = G AN7r (— m 2 .) 
rises with the quark mass because of the current-mass 
dependent pion pole location. This property is also vis- 
ible for <7 pTW in Fig. 5, albeit less pronounced, as the 
p— meson mass is non-zero in the chiral limit and there- 
fore varies over a much smaller range when evolving the 
current-quark mass. The shaded area in Fig. 6 indicates 
the threshold position Ma = M/v + m n where the decay 
channel closes, and the width of the band is again induced 
by the uj— dependence which enters mainly through the 
mass of the A, cf. Ref. [98]. 

We finally note that in determining the A — > Ntt tran- 
sition form factor we have neglected the pseudoscalar 
seagull terms which would appear in addition to the di- 
agrams displayed in Fig. 6. Judging from the smallness 
of the electromagnetic seagulls in the case of electromag- 
netic form factors, this approximation might be well jus- 
tified. The question of its validity can be settled by in- 
vestigating the A — > Ntt transition in the three-quark 
framework of Ref. [62] where all such missing contribu- 
tions, while no longer appearing explicitly, would be au- 
tomatically included. 



V. CONCLUSIONS 

We presented a calculation of the hadronic p — > tttt 
and A — » Ntt decays, as well as the pseudoscalar transi- 
tion form factor G AJV7r (G: 2 ), in the framework of Dyson- 
Schwinger and covariant bound-state equations. A con- 
sistent construction for the decay diagrams was imple- 
mented. The p — ¥ tttt transition was computed from 
the quark-antiquark Bethe-Salpeter equation in rainbow- 
ladder truncation whereas the A — ¥ Ntt transition was 
studied within the covariant quark-diquark model. All 
ingredients are determined self-consistently which leaves 
a phenomenological ansatz for the quark-gluon coupling 
as the only model input. 

The results in both cases compare well with experimen- 
tal and lattice data. The pirn coupling is underestimated 
by 15% and slowly rises with the current-quark mass, 
in agreement with lattice results. A similar observation 
holds for the A — > Ntt coupling which agrees also well 
with the experimental result. We find that G AJV7r (0) is 
practically independent of the current-quark mass. Con- 
sequently, the decay widths for p and A are mainly gov- 
erned by the available phase space. 

The present calculation provides a first step towards a 
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thorough investigation of hadron resonances and their de- 
cays within QCD. An important future direction in that 
respect would involve the implementation of explicit pinr 
and ANtt decay channels in the meson and A— baryon 
bound-state equations. Such an extension would repre- 
sent a considerable step beyond the rainbow-ladder trun- 
cation employed herein and contribute to a more realistic 
description of hadron resonances from their underlying 
dynamics in QCD. 
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and the polarization vectors of the p-meson are nor- 
malized to e^Ey = 6\y. In the case A — > Nir Eq. (12) 
entails 



M s , = u s A£ N7r <, = |jg. u s Q" <, , (A7) 



where Q is the pion momentum, and the nucleon and A 
spinors are normalized to 

^u s u s = 2M JV A+, 5>^ = 2Af A P^, (A8) 

s s 

respectively. This yields 

1 A/f A 

(A9) 



6M N 



with the spin sum 



<r = tr{A+(P / )A+(P i )} = l + Jl 



Mi 



(A10) 



Appendix A: Matrix elements and decay widths 



The strong decay widths in both cases are then given by 



For the decay of a particle with momentum p and mass 
M into two decay products with momenta pt and masses 
n%i, with p = YliPii the decay width is given by 



2M 



\M\\ 



(Al) 



where the squared transition matrix element |VW| 2 is av- 
eraged over the spins/polarizations, i.e. 



\M P ^\ 2 



\M 



(A2) 



and a sum over all final states is implicit. The phase- 
space factor / = k/{AitM) involves the quantity 



\/(M 2 - (mi + m 2 ) 2 ) (M 2 - (mi - m 2 ) 2 ) 



2M 



(A3) 



which in the rest frame of the decaying particle is given 
by k = \px | = |pa|; thus the decay width becomes 



r = 



\M\'< 



8ttM 2 

In the case p — > txtx one obtains specifically: 



(A4) 



M x = A^e» 



\-Mp-KTr\ 2 — g K p7r7r5p7T7r , (A5) 



where from Eq. (A3) and the kinematics in Fig. 3 one 
has 



m 2 /4 — m 2 = V P 2 



(A6) 



p _ K p7T7r 9p-rnv p _ <7 K AW7r 9&ntt / \ i i \ 



To extract the coupling constants from the matrix el- 
ements in Eqs. (10) and (16), one has to perform appro- 
priate momentum contractions and Dirac traces which 
yields: 



9k 2 



Appendix B: Color and flavor factors 



(A12) 



While the color and flavor traces in the quark 
DSE (2) and meson and baryon bound-state equations, 
Eqs. (3, 5, 7), have already been worked out in the main 
text, we still have to perform these traces for the p — > hit 
and A — > Ntt transition matrix elements (10) and (16). 

We work in the SU(2) /-isosymmetric limit and thus 
we have two degenerate quark flavors u — ( J ) and 
d = {i)- They transform according to the funda- 
mental representation of SU(2), whereas the anti-quarks 
u = ( J ) and d = ( i ) transform according to the com- 
plex conjugated fundamental representation. These can 
be used to construct representation matrices of mesons 
as quark-antiquark, diquarks as quark-quark and baryons 
as quark-diquark bound states. 

The 7r and the p-mesons discussed in this article are 
isovector states with three isospin projections, which can 
be labeled by the corresponding electric meson charge. 
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A possible representation is given by the matrices 



r + = \ud) = 



<J\ + ia-2 



ro 



\uu) — \dd) (T3 

ci - ic 2 
2 ' 



(Bl) 



where the cr, are the Pauli matrices, and the flavor ma- 
trices are normalized to tr{rj rj} = Sij. The color factors 
for the mesons are given by Sab, where A, B = 1,2,3 
denote the quark indices. 

The p°— meson in the triangle diagram of Eq. (10) can 
couple to the upper and the lower quark line. The color 
trace in both diagrams equals 3 whereas the Dirac traces 
yield an opposite sign: (A 1 ')^ 7r7r = -(A^)^, where A^ 
represents the expression in Eq. (10). The full color- 
flavor-Dirac trace of the p° — > -k + tt~ transition then 
yields 

~tr{r\_ rl_ r }3A T - tr{r^ r rl}3A l = 3\/2A T , 

i.e. the color-flavor trace of Eq. (10) is 3\/2. 

Similar to the mesonic case, for two degenerate flavors 
there exist four different diquarks, one of which is an 
isoscalar whereas the other three form an iso-triplet. The 
difference to mesons in flavor space is due to the different 
transformation properties of quarks and antiquarks. The 
representation for the diquark flavor matrices reads 

\u d) — \du) ioi 



so = 

Si = \uu) — 



V2 V2' 

1 +<T3 



S2 = 



\ud) + \du) o 1 

l-Q- 3 

~~2 ' 



(B2) 



s 3 = \dd) = 



The flavor factors for baryons in the quark-diquark pic- 
ture are given by the Clebsch-Gordan coefficients accord- 
ing to the respective diquark content of the baryon. For 
proton and neutron they read 



P = I u 

d 



d, - 
0, sj\d, 



u, 



(B3) 



where the first terms represent the isoscalar diquark con- 
tributions and the remaining three the contributions from 
the isovector channel in the same order as Eq. (B2). 
The A-baryons do not contain any contribution from the 



scalar diquark, and the corresponding Clebsch-Gordan 
construction yields 



A" 



u, 0, 



Id, 



u, 



A 1 



0, \Hd, 



A" 



(B4) 

0, 0, d\. 

Finally, the color factors are (sabc/V®) for the di- 
quark amplitudes and (Sac I V$) for the N and A quark- 
diquark amplitudes, where A, B are quark indices and C 
is the diquark index. 

In the case of the A ++ — s- pir + transition the three 
contributions in Eq. (17) yield the following flavor traces: 



3 ij 

£p]^WAt+. 



(B5) 



Only the axial- vector diquark contributes to the A (hence 
j = 1,2,3) whereas the nucleon has both scalar (i = 
0) and axial-vector diquark components (i — 1,2,3). 
Combined with the color traces +1 for the impulse- 
approximation diagrams and —1 for the exchange dia- 



grams the final result for the A +H 
matrix element in Eq. (16) reads 



pn + transition 



A = 



aA a 



A^A + ^3 A^A 



1 AAA 

2 Ji EX 



Vl ASA 
2 i4 -EX 



where the superscripts S and A refer to the scalar or 
axial-vector diquark content in the outgoing (left) and 
incoming (right) baryon amplitudes. 

The remaining processes are obtained accordingly by 
replacing p, A ++ and r + with the appropriate flavor 
structures of Eqs. (Bl) and (B3-B4). The bracket in 
the previous equation is identical in all cases whereas the 
prefactors become: 



• ±y I for the transitions A ++ — > p7r + , A — > rm ; 

• -§ for A+ pn°, A nvr°; 

• ±^ for A+ nn + and A -> pn~ . 

For every initial state one has to sum over all final states 
via Eq. (Al), hence one gets for both A + and A 



\ 



2 


2 


V2 


3 


+ 


3 



3 " 



(B6) 



and thus all flavor factors in the A7V7T system are the 
same. 
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